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Abstract 

In general relativity, it has been shown that the effective gravitational stress-energy tensor for 
short-wavelength metric perturbations acts just like that for a radiation fluid, and thus, in particular, 
cannot provide any effects that mimic dark energy. However, it is far from obvious if this property of 
the effective gravitational stress-energy tensor is a specific nature held only in the Einstein gravity, 
or holds also in other theories of gravity. In particular, when considering modified gravity theories 
that involve higher order derivative terms, one may expect to have some non-negligible effects arising 
from higher order derivatives of short-wavelength perturbations. In this paper, we argue this is not 
the case at least in the cosmological context. We show that when the background, or coarse-grained 
metric averaged over several wavelengths has FLRW symmetry, the effective gravitational stress- 
energy tensor for metric perturbations of a cosmological model in a simple class of f(R) gravity 
theories, as well as that obtained in the corresponding scalar-tensor theory, takes a similar form to 
that in general relativity and is in fact traceless, hence acting again like a radiation fluid. 

1 Introduction 

Our observable universe appears to be homogeneous and isotropic on large scales, but highly in- 
homogeneous on small scales. It is therefore considerably interesting to consider whether the local 
inhomogeneities can have any effects on the global dynamics of our universe, in particular, any effect 
that corresponds to a positive cosmological constant or dark energy. A number of authors have ex- 
plored this possibility of explaining the present cosmic accelerating expansion by some backreaction 
effects of the local inhomogeneities [THE]. Such a backreaction effect may be described in terms of 
an effective stress-energy tensor arising from metric as well as matter perturbations. 

In general relativity, a consistent expansion scheme for short-wavelength perturbations and the 
corresponding effective stress-energy tensor were largely developed by Isaacson [191120] . in which 
the small parameter, say e, corresponds to the amplitude and at the same time the wavelength of 
perturbations. Isaacson's expansion scheme is called the high frequency limit or the short-wavelength 
approximation. In this expansion, the dominant order of the Einstein equation with respect to 
this parameter e corresponds to the equations of motion for linearized gravitational waves in the 
ordinary perturbation theory, and is in fact divergent as e _1 . The next order of the expansion of 
the Einstein equation provides the Einstein equation for the background metric with an effective 
stress-energy tensor, which is essentially given as minus the second-order Einstein tensor averaged 
over a spacetime region of several wavelengths of metric perturbations. Since taking a derivative of 
perturbations corresponds, roughly speaking, to multiplying the inverse of the smallness parameter 
(or the inverse of the wavelength of perturbations), the effective stress-energy tensor consisting of the 
square of derivatives of the first-order metric perturbations can have large effects on the background 
dynamics. Furthermore it can be shown that the effective stress-energy tensor thus constructed 
is gauge-invariant, hence has a physical meaning. If the effective stress-energy tensor had a term 
proportional to the background spacetime metric, then it would correspond to adding a cosmological 
constant to the effective Einstein equations for the background metric, thereby explaining possible 
origin of dark energy from local inhomogeneities. It has been shown, however, that this effective 
gravitational stress-energy tensor is traceless and satisfies the weak energy condition, i.e. acts like 
radiation )2ip22j. and thus cannot provide any effects that imitate dark energy in general relativity. 

However, it is far from obvious if this traceless property of the effective gravitational stress-energy 
tensor is a nature specific only to the Einstein gravity or is rather a generic property that can hold 
also in other types of gravity theories. The purpose of this paper is to address this question in a 
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simple, concrete model in the cosmological context. Among many, one of the simplest of modified 
theories so far proposed is the so called f(R) theory, whose action is a generalization of the Einstein- 
Hilbert action to an arbitrary function, f(R), of the scalar curvature R. Since f(R) gravity contains 
higher order derivative terms, one can anticipate the effective gravitational stress-energy tensor to be 
generally modified in the high frequency limit. 

Over the past decade, cosmological implications of f(R) gravity theories have extensively been 
studied especially in the quest of finding an alternative cosmology to A-CDM model. A viable 
class of the f{R) theories is summarized in [23H3T] . Although it is desirable to examine all these 
cosmologically favored models, in the present paper, we will restrict our attention to the simplest 
model f(R) = R + cR 2 with an eye to applications to analyses of more generic cases, which are left 
for our future study. This model itself is not considered as a cosmologically favored modified gravity 
theory for describing the present accelerating universe, but has rather been introduced as a prototype 
of an inflationary universe model by Starobinsky [32]. However, this simple model can be viewed as 
the leading term truncation of a more generic class of f(R) theories that take an analytic form with 
respect to R around the vacuum solution R = and therefore provides, as the first step toward this 
line of research, a good starting point of our analysis. It would also be interesting to check whether or 
not a once-excluded model can possibly revive as a cosmologically favored model, due to the inclusion 
of the backreaction effects of local inhomogeneities. 

It is well-known that f(R) gravity is equivalent to a scalar-tensor theory, which contains the 
coupling of the scalar curvature R to a scalar field in a certain way 30 , 33 35 . The Brans-Dicke 
theory [36] is one of the simplest examples. Therefore our analysis can be performed, in principle, 
either (i) by first translating a given f(R) theory into the corresponding scalar-tensor theory and 
then inspecting the stress-energy tensor for the scalar field <j>, or (ii) by directly dealing with metric 
perturbations of the f(R) theory. One may expect that the former approach is much easier than the 
latter metric approach, as one has to deal with metric perturbations of complicated combinations of 
the curvature tensors in the latter case. Nevertheless we will take the both approaches. In fact, in 
the metric approach, by directly taking up perturbations of the scalar curvature R, the Ricci tensor 
R a b and the Riemann tensor R a bcd involved in a given f(R) theory, we can learn how to generalize 
our present analysis of a specific class of f(R) gravity to analyses of other, different types of modified 
gravity theories that cannot even be translated into a scalar-tensor theory, such as the Gauss-Bonnet 
gravity. 

The theory we consider in the present paper contains higher order derivative terms. As in the case 
of general relativity, we can consider short-wavelength metric perturbations with a small parameter 
e and expand the field equations with respect to e. In contrast to the Einstein gravity, the dominant 
part of the field equations for this theory is of order 0(e -3 ). In 0(e _1 ), we have equations of motion 
for linearized gravitational waves. In O (1) we obtain equations for the background metric with 
a source term arising from short-wavelength perturbations. This source term contains a number 
of higher order derivatives of metric perturbations. However, as one cannot have any meaningful 
notion of stress-energy for gravitational waves in a local sense (at least within a wavelength), we 
have to take a suitable spacetime average over several wavelengths. Also, since we are interested in 
backreaction effects on the cosmological dynamics, we assume that our background metric takes the 
form of the Friedmann-Lemaitre-Robertson- Walker (FLRW) metric. We also impose that in the limit 
to the Einstein gravity, i.e., f(R) — > R (as c — > 0), the field equations for our f(R) theory reduce 
to those for the Einstein gravity in the corresponding order of the expansion parameter. At this 
point, a number of terms that involve higher order derivatives of metric perturbations vanish by the 
spacetime averaging procedure and the assumption of the background FLRW symmetry. Eventually, 
besides the terms corresponding to the Isaacson's formula in the Einstein gravity, only a few terms 
that contain higher order derivatives of metric perturbations can remain in the effective stress-energy 
tensor for short- wavelength perturbations in our modified gravity. Furthermore, the resultant effective 
stress-energy tensor is shown to be traceless as in the Einstein gravity case. 

We briefly comment on the previous work along the similar line. Since the effective stress-energy 
tensor can be used to measure energy flux carried out by gravitational radiation from astrophysical 
sources, such as inspiral binary systems, it can be used to test various modified gravity theories by 
the near future gravitational wave detectors [37H45] . For this purpose, the effective stress-energy 
tensor for gravitational radiation has been derived by Sopuerta and Yunes [46] by applying Isaacson's 
scheme to the field equations of dynamical Chern-Simons theory. A more general formalism to 
compute the effective stress-energy tensor from the effective action, which can apply to a wide class 
of modified gravity theories, has been proposed by Stein and Yunes [47]. As a concrete example, 
the formula has been applied to dynamical Chern-Simons gravity as well as theories with dynamical 
scalar fields coupled to higher-order curvature invariants. It has been shown that in these modified 
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theories the stress-energy tensor for gravitational radiations reduces, at future null infinity, to that 
in the Einstein gravity. Berry and Gair [48] also have derived the effective stress-energy tensor for 
gravitational waves in the f(R) gravity which is analytic around the vacuum R — 0. However, since 
the main concern in these studies is mainly to test alternative theories in the astrophysical context 
by using gravitational wave detectors, the formulas mentioned just above have been formulated for 
asymptotically flat spacetimes (as energy flux of gravitational waves needs to be evaluated at future 
null infinity) and therefore do not appear to apply to cosmological models. In contrast, our analysis 
will proceed by exploiting the cosmological setup that our background metric possesses the FLRW 
symmetry. 

In the next section, before going into the effective stress-energy tensor in modified gravity theories, 
we will first briefly summarize the high frequency limit in general relativity. In Sec. [3] we consider 
the high frequency limit in f(R) gravity theory. Based on the Isaacson's scheme we expand the 
field equations for f(R) = R + cR 2 theory and first derive the general expression of the effective 
stress-energy tensor for gravitational perturbations in our f(R) gravity. Then, assuming that our 
background metric has the FLRW symmetry and also that the resulting equations reduce to the 
corresponding equations for the Einstein gravity in the limit c — » 0, we see that the effective stress- 
energy tensor whose expression is significantly simplified, is in fact traceless as in the Einstein gravity 
case. As briefly mentioned above, when a given f(R) gravity is translated into the corresponding 
scalar-tensor theory, the scalar field <j>, which expresses an extra-degree of freedom in the f(R) theory, 
possess a non-trivial potential term. [Compare with the earlier work by Lee [49] on a computation of 
the effective stress-energy tensor in a scalar-tensor theory with vanishing potential term.] In Sec. [4] 
we will make sure that the effective stress-energy tensor in Brans-Dicke theory is consistent with that 
in our f(R) gravity. Sec. [5] is devoted to summary and points to future research. 

Our signature convention for g a b is ( — ,+,+,+). We define the Riemann tensor by R a bc d ^d = 
2V[ a Vi,]o;c and the Ricci tensor by R a b = Racb c as in Wald's book [50] . 



2 High frequency limit in general relativity 

In this section we introduce our notation by recapitulating Isaacson's expansion scheme for short- 
wavelength gravitational perturbations in general relativity. 

Let g a b be the metric with linear perturbation h a b\ it is described by g a b = g™ + h a b with g™ 
being the background metric including the backreaction from perturbations. The amplitude of h a b 
is of order h a b ~ O (e) with e being the small parameter, which also corresponds to the wavelength 
A of perturbations compared with the background characteristic curvature radius, L. The order of 
derivatives of h a b are 

V ai V Q2 ■ ■ ■ V am h bc ~ O (jjjfiz) ~ O (e 1 -™) , (2.1) 

where V a denotes the covariant derivative with respect to gf b \ so that V a g£° = 0. We may bear in 
mind perturbations of the form h ~ esin(j;/A) and X/L ~ O (e). In what follows we normalize L ~ 1. 
The inverse metric takes the form g ab = g {0)ab - h ab + h a c h cb + ■■■, where h ab = g (0)ac g (0)bd h cd . 
There is the general relationship between the Ricci tensor of g a b and that of gf b \ namely 

Rab = Rab[g {0) ] + 2V [e C£j a + 2C c d[c ct ]a , (2.2) 

where C£ c = r£ c — r£ c [</ '] = g ad (Vbgdc + V c <?m — Vdgtc) /2. Since R a b contains terms such as those 
schematically expressed as g _1 Wg, Vg~ x Vg and g~ 1 g~ 1 \7g\7g, we can find 

R%>[h]~0{e n - 2 ) , (2.3) 
where n is the number of h ab included in Rab- We also find 

R™[h]~G™[h]~0(e n - 2 ) , (2.4) 

where R = g ab R a b, G a b = Rab — g a bR/2 is the Einstein tensor, and R^lh] and [h] do not contain 

Rabi9 l °>]. 

The Einstein equation is 

Gab = Rab - \gabR = k % T™ , (2.5) 
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where k 2 = 8ttG and is the stress-energy tensor for the background matter fields. In the following, 
for simplicity, we focus on metric perturbations and ignore perturbations of the matter fields in T^' . 
We can find the dominant terms, O of the Einstein equation as 

or simply 

R%[h] = 0. (2.6) 

This is equivalent to the equation for linearized gravitational waves in the ordinary perturbation 
theory. Next, in the order of O (1), we find 

G ab [g {0) ] = * 2 Ti°' + K 2 T a f , (2.7) 
where the effective gravitational stress-energy tensor, , is given by 

« - -(<?£> [ft]) 

= (iv a ft™V 6 ^ T ^ . (2.8) 

Here and in the following (■ • •) denotes taking a spacetime average over several wavelengths of per- 
turbations. Here the indices are raised and lowered with g^ ah and gS ■ For the expression of the 
third-line, the transverse-traceless gauge, V c h c a = = h c c , denoted by h TT , and (|A.6|I have been 
used (see also Appendix. fSJ. We can check that is traceless, i.e., it acts like radiation: 

n 2 T eSa a = 0, (2.9) 

from (|2.6[) and (|A.2[) , Thus, in particular, it cannot provide any effects that mimic dark energy in 
general relativity. For more mathematically rigorous treatments of short-wavelength perturbations 
and the effective stress-energy tensor, see [211122] , 

The effective stress-energy tensor (|2.8[) can be shown to be gauge-invariant [21]. In fact, the 
expression of the right-hand side of (|2.8|) is given by manifestly gauge-invariant part of hj^ . For 
this purpose, one can introduce the polarization tensors eil' , as usual, and decompose the metric 
perturbation accordingly hj^ = e^h^ + e^ft' x \ In the cosmological context, one is concerned 
with the Friedmann-Lemaitre-Robertson- Walker (FLRW) metric, 

ds 2 = -dt 2 + a(t) 2 j id dx i dx j (2.10) 

with da 2 = ~iijdx % dxi being the metric of 3-dimensional constant curvature space. So, it may be more 
convenient to impose the transverse-traceless condition with respect to this FLRW time-slicing, i.e., 
^oo T = 0, h^J = 0, V a h TTa i = 0, h TTl i = 0. This condition completely fix the gauge freedom and 
(12. 8p is written by the gauge invariant variable hj^ as ( jVah TTlJ : V(,/iJ T ). 



3 High frequency limit in f(R) gravity 

The general action for the f{R) gravity is given by 

where Cm is the Lagrangian for matter fields, such as perfect fluid in the cosmological context. 
Varying this action with respect to the metric, we have the field equations 

GH H) = Gab + FR ab - ^g ab f - VaV b F + g ab g cd V c V d F = k 2 T^ , (3.2) 

where f = f — R, F = df /dR, and denotes the matter stress-energy tensor. 

The field equations in the f{R) gravity have terms consisting of higher order derivatives of R, and 
the order of those derivatives are higher than that of R: 

V ai Va 2 ■ ■ ■ Va m R {n) [h] ~ O ( e n ~ 2 - m ) . (3.3) 
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Therefore it is expected that the effect of the short-wavelength approximation would be enhanced. In 
order to see whether this is the case, from now on we restrict our attention to the following concrete 
model 

f(R) = R + cR 2 , (3.4) 
where c is a constant. This model has been considered for the first time in the context of inflationary 
universe [32]. The field equations are 

G f a [ R) = G ab + 2c (RR ab - ±g ab R 2 - X7 a V b R + g ab g cd V c V d R^j = k 2 t£> . (3.5) 

As in Isaacson's formula reviewed in the previous section, we expand the above equations with 
pect to the small par 
the following equations 



respect to the small parameter e. Then, the dominant part is of the order O (e 3 ) , in which we have 



V a V 6 i? (1) [h] - g^UR™ [h] = . (3.6) 
By contracting with g^ ab , we immediately have 

V a V b R w [h] = 0. (3.7) 

Next, for the order O (e 2 ), we have 

R W lh]R$[h] - \g™ ( R V[h]) a - V a V b R™[h]+g$nR™[h] = 0. (3.8) 

By dotting with g (0)a6 , we have 

DR (2) [h] = 0, (3.9) 

V a V b R^[h] = R w [h]R^[h]-\g^(R w [h]) 2 . (3.10) 

Note that since we are working in the short- wavelength approximation, we find g^ ab R 1 ^ [h] = i?' 1 ' [h] 
in 0(e _1 ), which is different from calculation in ordinary perturbation theory, where in general 
9 {0)ab Ra b [h] £ R w [h]. For O (e" 1 ), we have 

R$ [h] ~ ~9%R W [h] + 2c { (R[g^] + R^ [h]) R$ [h] + R™ [h] ( Rab [ g ^] + [h]) } 

-f {Rl9 m ] + R {2) [h]) R w [h] + h ab (#«M) 2 } 



-2c [V a V b R (3) [h] - Rcd[g ( °>]V a V b h cd ) 

+2c{g^ (aR^[h]-R cd [g m ]ah cd ) - g^ h cd V C V d R {2) [h]} = 0. (3.11) 

Again, by dotting with g^°' ab , we have 

R w [h] + ^h a a (i? (1) [ft]) 2 - 6c (aR (3) [h] - R cd [g (a) ]Uh cA j + 8ch ab V a V b R {2) [h] = 0. (3.12) 

In order O (1), as in Isaacson's formula in general relativity, we have the field equations for the 
background metric with the backreaction source term: G{£ '[</ '] = K 2 T^°' + K 2 T^f , where 

« = - \ (g<£RV[h] + h ab RV[h]) 

+2ci [ R^[h]R^[h]+R[g^]Ri 2 b ) [h] 

+R™[h] (iU[s (0) ] + (fl< 3 >[fc] - h cd RM 0) }) R$[h] 

~\ \9 a l } {2R W lh}R^[h] - 2R^[hW d RM 0) ]+2R[9 i0) ]R {2) [h] + (^[fc])''} 

+2h ab (R{g W ]+R (2> [h])R W [h] 
+2c{ (-si°> cd + W 0)cd ) (y c V d R {3) [h] - i? e/ [3 (0) ]V c V d ft e/ ) 

+ (9abh c eh ed - h ab h cd ) V c V d i? (2) [/i] j V (3.13) 
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This is the expression of the stress-energy tensor for short-wavelength metric perturbations on the 
generic background metric g a °^ in our f(R) gravity. 

From now on, we consider in the cosmological context. We assume that our background is spa- 
tially homogeneous and isotropic, that is, our background metric possesses the FLRW symmetry and 
therefore takes the form of (|2.10[) . Then, thanks to this background symmetry we can explicitly 
solve equations of the form X7 a V b S(t,x) — 0, such as (|3.7|) (see Appendix. [Bj. Equation (|3.7|) (the 
equations of motion of O (e -3 )) is solved to yield 

R w [h]= const. (3.14) 

Taking the average, we find 

R w [h] = const. = (const.) = ^R m [h]) =0. (3.15) 
Then, the equations (|3.10[) (those of O (e~ 2 )) become 

V a V b R (2) [h] =0. (3.16) 
Again using the result in Appendix. [B] we find 

R {2) [h] = Si = const . (3.17) 
By using (|3.15|) and (|3.16[) . the equation (|3.12[) (of O (e - )) immediately yields 

□i? (3) [h] - R ab [g m }ah ab = , (3.18) 
and the equations (|3.11l) (of O (e -1 )) reduce to 

(l + 2cR[g (0) ] + 2cSl) R% [h] = 2c (v a V b i? (3) [h] - R cd [g (0) ]V a V„/i cd ) . (3.19) 
The effective stress-energy tensor (|3.13p is then expressed as 

« = J(l + 2cR[ g ^])R^[h]- 1 -g^Si 

+2c ' 



{ {R^ [h] - h^R cd W 0) \) R$ [h] + Si (W>] + R%>) } 
- C 1 9i° b ) {2R[g (0) ] + s 1 )s 1 

+2c (-<?£° 6 V d ) (v c V dJ R (3) [fe] - R af [gW]V c V d h ef )\ . (3.20) 

Since the effective stress-energy tensor in the R 2 model should reduce to that in general relativity 
when c = 0, we choose Si (= R {2) [h\) to be 0. Then, ((^20|) becomes 

« = -/(l + 2cRfo<°>]) R^[h] + 2c (R&[h] - h cd RM 0) \) R$[h] 
-g$ (l + 2cR\ g ^])h^[h]\ 

= -^(l + 2cR[g^]) R$[h]+2c (RW[h] - h cd R cd [g (0) ]j R$[h]j , (3.21) 

where we have used (13 . 19 P in the first equality above, and 

[h*R$[h])=0 (3.22) 



in the second equality so as to make the above expression compatible with that of general relativity 

in the c = case. The expression, (|3.21[) . is our main result of this section. From R^[h] = and 
(f3T22]> . we see 

(g mab Rfhh]) = (RW[h] + h« b R$[h]) 

= . (3.23) 
Then using this and R (1> [h] = 0, we can find that K 2 T„f is in fact traceless: 

K 2 T eSa a =0. (3.24) 



G 



4 The high frequency limit in scalar-tensor theory 



In the previous section, the scalar curvature R and the Ricci tensor R a b are taken up directly in the 
metric formalism of the f(R) gravity. It is well-known that any f(R) gravity theory is included in 
Brans-Dicke theory, which is one of the simplest examples of scalar-tensor theory 30 ,33 35 . In this 
section, we will see that the results obtained in the previous section are indeed consistent with those 
obtained within the corresponding scalar-tensor theory. 
The action of Brans-Dicke theory [36] is 



V a 4>V a <t>-V{cj>) \ + J d 4 x£.i 



where ujbd is a constant called the Brans-Dicke parameter and <f> is a dimensionless scalar field, and 
Cm denotes the Lagrangian for matter fields, which can in general couple to the metric g a b as well as 
the scalar field cj>- Then the equations of motion for <j> and g a b are, respectively, obtained as 

<t>G ab - ^ (V„^V„0 - ?f V c 0V c <^) - V a V b <t> + g ab (V d V c V d + V(<ft) = k 2 t£> , (4.3) 

where T^ ^ and T^°' are the stress-energy tensor for matter fields obtained by taking variations of <f> 
and g a b, respectively. 

The f(R) gravity of the metric formalism, (|3.1[) . can be cast into the form of the above Brans-Dicke 
theory by setting 

, „ (m _df(R) F(R)R — f(R) 

<t> = F(R) = , ojbd=0, V = . (4.4) 

In this case, as one can find R — 2d^,V = 0, the equations of motion for <f> and g a b just given above 
become respectively 



1 

2^(0) 



□0__L V >V^ = K 2 T^ 0) ab , (4.5) 



G^b = <t> [Rab ~ -gabRJ - V a Vb<P + Qab V c V d + V (<f>) ) = K^T^ . (4.6) 

From now we consider short-wavelength perturbations for <j>: <j> — 4>o + Scf>. We also assume that 
there is no coupling of matter fields with the second-order derivatives of <f>, so that there are no 
non-vanishing terms of order 0(e _1 ) in the stress-energy tensor for matter fields. Then, the equation 
of motion for of O (e~ ) is 

as</> = , (4.7) 
and the equations of motion for g a b of O (e _1 ) arc 

<h (r$ [h] - \g^R {1) [h]j = V„ V b 5<t> - g$>n54> . (4.8) 

Contracting with g^ ab , we have 

R w [h] = -^-050 = 0, (4.9) 

00 

where we have used (14.71) . From this equation, we can immediately find 

^VaVfcty. (4.10) 

The equations of motion in O (I) are given by G^J [g {0) , </> ] = K 2 T^ } + k 2 T^ , where 

k 2 T c J = - (foR.™ [h] + 8<t>R)2 [h] - g^h cd V c V d 5<p) . (4.11) 

Here we would like to emphasise that so far we have made no assumptions concerning the form of 
f(R) or the symmetry of our background metric g a ° b , the above expression, (|4.11|) . applies to the 
generic f(R) theory with an arbitrary background metric. 

If we restrict the form of f(R) to be (|3.4[) . then by inspecting the expansions <j> = <f> + Scf> + ■ ■ ■ 
and F(R) = 1 + 2cR = (l + 2cR[g (0) ]j + 2c - h cd R cd [g (0) ]j + ■ • ■ , we find 

00 = 1 + 2cR[g m ] , Sep = 2c (# (3) [h] - h cd R cd [g m ]) . (4.12) 



7 



Using these and (|4.10[) . we have 



**Tg = - ((l + 2 CJ R[ 3 (0) ]) R^[h] + 2c (R^[h] - h cd RM° ] ]) B$[h] ~ 9% ^ h cd R^[h]) .(4.13) 

Now we work on the cosmological situation so that the background metric has the FLRW sym- 
metry. Provided that the limit c — > should reproduce results in the case of the Einstein gravity, we 
finally obtain 

« = - ( (l + 2cR[ff<°>]) [h] + 2c (fl< 8 > [h] - h cd RM 0) \) R$ [/»]) , (4-14) 

where we have used (|3.22|) . derived under the FLRW symmetry in Sec. [3] We see that the expression 
(|4.14[) above is precisely the same as (|3.21[) derived within the metric formalism of the f(R) gravity. 
This verifies our methods of Sec.[3]for dealing with short-wavelength perturbations of the f(R) gravity 
within the metric formalism. 

5 Summary 

We have addressed the effective gravitational stress-energy tensor for short-wavelength perturbations 
in the simple class of f(R) gravity of R 2 type in the cosmological context. As in the Isaacson's 
formula for the Einstein gravity reviewed in Sec. O we have obtained the field equations for the 
background metric with a backreaction source term T£§ in order O(l) of the small parameter e. 
Reflecting the fact that our f(R) theory contains higher order derivative terms, the source term or 
the effective stress-energy tensor takes, as given in (|3.13[) , quit a complex form that contains, in 
principle, terms of fourth order derivatives, schematically expressed as (V/iV/iV/iV/i). The resultant 
expression, (|3.13[) . of the effective stress-energy tensor, in fact, applies to any background metric 
9ab' Until this point, no symmetry assumption on the background metric has been used. Then, by 
imposing that our background has the FLRW symmetry, we have derived our effective stress-energy 
tensor for short-wavelength metric perturbations in cosmological models. At this point, thanks to the 
background FLRW symmetry and the spacetime averaging over several wavelengths, the expression 
of our effective stress-energy tensor has been significantly reduced to have the simple form, (|3.21l) . 
We have also shown that the obtained effective stress-energy tensor is traceless, so that it acts like a 
radiation fluid as in the Einstein gravity case and thus, in particular, cannot mimic dark energy. 

Since any f(R) gravity theory is known to be equivalent to a scalar-tensor theory, we have cast 
our f(R) theory into the corresponding scalar-tensor theory. Then, within the scalar-tensor theory, 
we have derived the effective stress-energy tensor for short-wavelength perturbations of the scalar 
field and checked consistency with the stress-energy tensor obtained within the metric perturbations 
of the original f(R) theory. 

Although we have focused on the R 2 model especially about the FLRW background, in Sec. [4]we 
have pushed forward our calculations with a general f(R) gravity about an arbitrary background as 
far as possible, and have not used the property of the R 2 model about the FLRW background, up to 
(14. lip . We can immediately note that (|4.11[) does not involve any terms of fourth order derivatives 
but has only terms of the square of first order derivatives of perturbations h a b and 8(f). This result 
obtained within the framework of the scalar-tensor theory indicates that the higher order derivatives 
could vanish also in the metric framework of general f(R) gravity theory for a generic background. 
However, to see whether this is indeed the case needs further involved calculation, and is beyond the 
scope of this paper. This is left open for future study. 

Our formulas derived in Sec. [3]deal directly with the scalar curvature R and the Ricci tensor R a b, 
and therefore should be able to apply to similar analyses of other modified gravity theories which 
contain higher order curvature terms composed of R, R a b, and R a bcd and which cannot even be cast 
in the form of a scalar-tensor theory. It would be interesting to consider an extension of our present 
work to a wide class of modified gravity theories with high-rank curvatures. 
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A Perturbation formulas in Transverse- Traceless gauge 



The conditions of transverse-traceless waves are X7 a h a b = h a , 
axe 

cT[h TT ] = 0, (A.l) 

= -\^h T J (A.2) 
R%[h^] = 2V lc C^[h^]-C^[h^]C^[h^] 

= --jh Vc\Vbh da +V a h bd ~V d h ba J 

+i v b [h TTcd (v c hJJ + v a h T J - v d h T j) } 

-\{v d h TT \ + V^ TT / - V c ftJ 6 T ) (V c h™ a + V a h TT / - V d h??) , (A.3) 

and the scalar curvatures for h^ b are 

R w [h TT ] = g^R^h™] 

= 0, (A.4) 

[h TT ] = g {0)ab R a 2 b ] [h TT ] - [h TT ] 

= ^V a h™^X7 -^ a h TTbc V c h T J + h TTab Dh T J , (A.5) 



where we have used [V a ,V b ]h]:J = R c \ b [g {0) ]hJJ + R d e ab \g {0) }h^ = 0(e). From flU]) and l|A3 
we find Dft^ = 0. Using these we find 



i?f> TT ]) = -iv a ft TT «V 6 ^ d T , (A.6) 
J R (2) [/i TT ]\ = 0. (A.7) 



B Solution of V a V b S{t, x ) = 

We solve the equation 

V VbS(t,2) = 0, (B.f) 
such as (|3.7|) . in the FLRW background spacetime: 

ds 2 = -dt 2 + a 2 (t)ji d dx i dx j . (B.2) 

Since we are interested in the expanding universe, in what follows we assume that the scale factor 
a(t) is dynamical, i.e., a(t) 7^ const. The (0,0), (0, i) and components of (|B.f I) are 

S(t,x) = 0, (B.3) 

diS(t,S)-^-diS(t,S) = 0, (B.4) 
a(t) 

did j S(t,x)-a(t)d(t)y l] S(t,x)-r k l ^ ) d k S = 0, (B.5) 



where the dot denotes the derivative with respect to the cosmic time. The solution of (|B.3|I is 

S(t,x) = Cl (x)t + c 2 {x) , (B.6) 
where c\(x) and C2(x) are arbitrary functions of x. Since Equation (|B.4|) becomes 

^-t)d i c 1 (x)=9M£)- (B.7) 
a(t) J 

Therefore, either 

ci,C2 = const. (B.8) 

or diC2/diC\ — const.. In the latter case, by shifting t — ¥ t — const., we have a/a = 1/t, and therefore 
have a(t) oc the behavior of the background FLRW universe is determined. This is, in our present 
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context, too restrictive, and for this reason, we should take the former case; ci,C2 = const.. Then, 
equation (|B.5|) for i — j becomes 

a(t)a(t)j i:j ci = , (B.9) 

which immediately implies 

ci=0. (B.10) 

Therefore we find the solution of (IB. II) to be 



S = const. (B.ll) 
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